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SYSTEM AND METHOD FOR GEOREFERENCING MAPS 

CROSS REFERENCE TO RELATED APPLICATIONS 

This application is a continuation-in-part ("CIP") of US patent application, Serial 
5 Number 09/537,849, filed 3/29/00 and shares specification and- figures with commonly.; 
assigned, copending US Patent Applications entitled "System ahd'Methodfor Performing " 
Flood Zone Certifications" (Serial Number 09/537,161 filed March 29, 2000); a similarly 
titled CIP of the '161 application, Serial Number , filed ; "System and 

Method for Synchronizing Raster And Vector Map Images" (Serial Number 09/537,162 
10 filed March 29, 2000); and a similarly titled CIP of the '162 application, Serial Number 
, filed ; all of the foregoing being incorporated by reference hereinto. 

FIELD OF THE INVENTION 

The present invention generally relates to geographic map processing and, in 
15 particular, to a system and method for georeferencing originally un-georeferenced 
paper maps, that is, for assigning latitude and longitude references to paper maps 
originally lacking them (or any other geographic referencing). More particularly, the 

present invention relates to scanning un-georeferenced paper maps to produce un- 
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georeferenced digital raster maps, and then associating the digital raster maps with a 
geographic coordinate system. The present invention can be advantageously used to 
georeference flood maps published by FEMA. 

BACKGROUND OF THE INVENTION 

5 FEMA FLOOD MAPS 

The Federal Emergency Management Agency (FEMA) publishes a large set of 
"flood maps." These paper maps are the official, Federally sanctioned source for 
determining flood risk. The flood maps delineate regions that are assigned a flood zone 
designation from among a limited set of possible designations. 

10 For example, a region designated "Zone A" on a flood map is a high risk area, 

which FEMA believes has at least a 1 in 30 chance of flooding during any given year. A 
region designated "Zone V" combines the properties of a "Zone A" with the additional 
risk of wave action or flowing water. In the foregoing context, performing a "flood zone 
determination" (also known as a "flood hazard determination" or a "flood certification") 

15 relative to a location entails a two-step process: (1) determining, based on FEMA flood 
maps and other materials, the FEMA-designated flood zone for the location in question, 
and, then, (2) completing required paperwork to report this determination (along with 
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certain other required information) to interested parties, such as home buyers, their 
lenders, and real estate brokers. 

FEMA flood maps come in a variety of styles, but broadly speaking they can be 
divided into two types of maps, namely, index panels ("IP"), or index panel maps, 
5 Figure 5, and regular panels Q X RP"), or regular panel maps, Figures 1-4 

Referring to Figure 1, a regular panel map 10 shows a relatively small geographic 
area at a fairly high level of detail. Each regular panel 10 depicts various "features" of 
the relevant area, these features include (without limitation) streets 12,13 railroad 
tracks 14, streambeds 16, lakes, intersections 18,20 of the foregoing, and detailed 
10 delineations of the community borders and the boundaries of various extant flood zones 
(not shown in Figure 1, but see Figure 2) and their classifications. 

Figure 2 depicts a regular panel map 10 in which flood zones 22 are depicted. 
As is typical, the flood zones 22 are denoted or "flagged" as variously shaded regions 
contained within irregular curvilinear boundaries (the "flood zone boundary")- Labels 
15 30 on regular panels 10 designate the types of flood zones that have been shaded. 
With the exception of regular panels 10 described as "not printed" ("PNP") and panels 
10 presenting flood zone-designated communities, described in greater detail below, all 
flood zone determinations are made by reference to appropriate regular panel maps 10. 
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Referring to Figures 1 and 2, each regular panel map 10 is bounded by a 
"neatline," that is, a rectangular frame or border 32 surrounding the geographic area of 
interest. The geographic area of interest may not completely "fill" the neatline 32, as 
indicated at 34 in Figure 1. For example, the actual geographic area contained within 
5 the neatline 32 of a regular panel 10 might consist of (a) a portion 33 of a community 
and the designated flood classification of the portion and (b) unincorporated or other 
portions 34 which surround, are surrounded by, or abut the community, arid which 
have no flood classification depicted on the panel 10. 

Referring to Figure 3, if FEMA has designated a regular panel maplO as being 

10 applicable to determining the flood zone classification of only an incorporated 
community, designated at 42 in Figure 3, then the regular panel 10 cannot be used to 
determine the flood classification of any depicted unincorporated areas, designated at 
34,36,38. According to FEMA protocol, the foregoing obtains whether or not flood 
zones for the unincorporated area are shown on the regular panel 10. Indeed, a FEMA 

15 panel map 10 often does not show flood zones for areas outside of the designated 
region for the map 10 (the designated region 42 is a portion of the community in this 
example). Sometimes the area outside of the designated area is just shown as a blank 
area on the panel, as at 34,36,38. Regardless of how much or how little detail is shown, 
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by convention, a flood determination can be made using a panel only in the region 
designated by FEMA. 

The "inset" 44 of a panel 10 means any geographic region that is both shown on 
the panel and as to which FEMA mandates the flood zone designations set forth on the 
5 panel to be applied. A panel inset 10 may be self-connected (connected with all parts 
of itself) as in Figures 1-3, or, as shown in Figure 4, it may be split into, two or more 
self-connected regions 46,48 by intervening non-designated areas 50 on the panel -10 
and may include other non-designated areas, such as those shown at 52 and 54. The 
"insets" of a panel map 10 may comprise a collection of disjoint, self-connected regions 
10 46,48 which collectively contain all of the designated area of the panel 10, and 
themselves constitute a subset of the geographic region contained within the neatlines 
32. 

Thus, from the foregoing, an actual map panel 10 consists, in general, of three 
parts: (1) the inset(s), (2) other (non-designated) geographic regions within the 
15 neatlines 32, but not included within the inset(s), and (3) areas of the panel which are 
outside of the neatlines 32, and generally show various types of textual information 
such as a map legend, the map panel number, and panel effective date, etc. It is true 
that, in the ideal case, the insets 46 and 48 should include only designated regions of 
the panel 10, but as a practical matter it is sometimes convenient for FEMA to include 
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some non-designated areas within the insets, so as to simplify the process of gathering 
inset border data. 

Referring to Figure 5, an index panel map 60 shows a considerably broader area 
than a regular panel 10, but at a lesser degree of detail, and usually without the flood 
5 zones being indicated thereon. Typically, an index panel 60 shows a broad geographic 
„ view of the locations of the various regular panels 10 that make up an entire 
community, or large portion thereof. It should be noted that;^hile FEMA ^defines" 
communities, these FEMA-defined communities often do not necessarily correspond to 
one's intuitive sense of what a community is. 

10 The index panel 60, generally depicted in Figure 5, provides its broader, less 

detailed overview by displaying the neatline boundaries 32 of each of its regular panels 
10 against a background of major geographic features (e.g. highways, lakes, etc.). 
Each of the regular panels' 10 neatline boundaries 32 is labeled with the corresponding 
number of its regular map panel, as represented by W XX," W YY" and n ZZ" in Figure 4.. 

15 Index maps 60, in effect, show how to "put together" the panel map 10 "jigsaw", that 
is, how to associate panel maps 10 with each other to form a grid thereof, grid lines 62 
thereof being comprised of the aligned and/or overlapping neatlines 32 of adjacent 
panels 10. 



6 



108344-00010 




Until the advent of the present invention, the primary purpose of the index maps 
60 has been their use as an aid in manually determining which of the regular maps 10 
should be consulted in order to make a flood zone determination. For example, if a 
property is known to be located in Community A, then a map analyst might first locate 
5 the property on the Community A index map 60. The analyst could immediately 
determine --by examining the neatline borders 32,62 and the labels W XX," W YY," W ZZ" 
displayed on the index map 60— which regular panel map 10 to consult. Since there 
are hundreds of panel maps 10 in the index map 60 of various communities, 
identification of a particular panel map 10 by examining an index map 60 can be 

10 effected more quickly than serially examining large numbers of panel maps 10. The 
present invention provides other, much more efficient ways of determining the correct 
regular panel map 10 applicable to a particular location. 

Index maps 60 also contain neatline borders 32 for panel maps 10 that are not 
printed as such, and which are designated ''panel not printed" or "PNP." As this name 

15 suggests, no actual panel map 10 is ever printed for a PNP. Instead, the relevant index 
map 60 displays the flood zone classification for the region contained within the PNP 
neatline borders 32. Accordingly, In the case of a PNP, the flood zone classification is 
determined from the index map 60, rather than from a regular panel 10. 



108344-00010 




In addition to the foregoing techniques, flood zone determinations can also be 
made on the basis of FEMA's denoting certain communities as "flood zone designated 
communities." In this event, FEMA mandates that a denoted community have the same 
flood zone classification throughout. As a consequence FEMA does not publish any 
5 flood panels —either regular 10 or index 60— for the community. Property located 
anywhere in thatcommunity assumes the,flood..zone classification of the community. 

DIGITAL RASTER MAPS 

A digital raster map image, Figure 2, may be thought of as a number of points 
selectively present on —or missing from— each X,Y coordinate of a Cartesian display. 

10 Points that are present are Ousually called "pixels." Besides having an address at an 
X,Y location, a pixel may have a "value," usually a numeric quantity designating a 
specified color or grey-scale value for the pixel. In this latter event, a "missing pixel" 
may be represented by a pixel having a value pursuant to which it is colored white. 
The X,Y location and value of each pixel is data that may be stored. A digital raster 

15 map may be created by optically, digitally scanning a conventional paper map and 
thereafter using a computer system to define and/or display it in terms of its pixel 
locations and color values. Commonly used formats for storing scanned images 
currently include TIFF, JPEG, and PNG. Any FEMA regular panel 10 or index panel 60 
may be represented as a digital raster map. 
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DIGITAL VECTOR MAPS 

A vector map, Fiogure 6, may be viewed as a "virtual" map that comprises, in 
effect, a collection of stored data. The stored data define the locations of plural nodes 
and plural straight line segments interconnecting the nodes. The stored data defining a 
5 vector map represents geographic information in the form of collections of points —the 
nodes— and line segments interconnecting the nodes to form polygons, i.e., closed 
figures comprised of line-interconnected nodes. These points, lines, and polygons may 
have additional information associated with them. For example, the stored data may 
indicatre that a particular point or node is the location of the capital building, or that a 

10 selected polygon describes the border of Lake Placid. The polygons are surrogates for, 
or representations of, the actual curvilinear configuration of various geographic 
features: streets, streams, lakes, cities, counties, etc. Typically, the features —nodes 
alone or two or more interconnected nodes— are presented in terms of 
latitude/longitude (or lat/lon) values. If this or another earth-referenced system is used 

15 to locate the nodes of the map, the map is said to be georeferenced. A visually 
sensible map may be "drawn" from the stored data 

In view of the fact that georeferencing of digital raster map images plays a key 
role in the present invention, the precise meaning of this term, for the purposes 

hereofand as used herein, is set forth below. 

9 

108344-00010 



Georeferencing 

A raster map image is said to be "georeferenced" if there exist mathematical 
functions, f and g, such that one function can be used to convert the X,Y coordinates of 
the pixels describing the map image to corresponding latitude/longitude coordinates 
and the other function can be used to convert the latitude/longitude coordinates of the 
pixels of the image to their corresponding X,Y coordinates. In other words, f and g 
effect the following: 

1. If (x, y) represents a pixel location on a digital raster map image, then 
f(*/ y) = (Lon, Lat) represents the longitude and latitude of the 
corresponding physical location; and 

2. If (Lon, Lat) represents a physical location that lies within the region 
covered by the map, then g(Lon, Lat) = (x, y) represents the point on the 
digital raster map image that corresponds to that longitude and latitude. 

In the previous paragraph, (x,y) defines a point pursuant to the natural internal 
coordinate system of the raster map image. Most digital vector map images use 
longitude and latitude as their internal coordinate system. Accordingly, most digital 
vetor maps can be considered to be trivially georeferenced per se. Most digital raster 
map images use the pixels of its image as a kind of natural coordinate matrix. Such a 
raster map image is not inherently trivially georeferenced and requires the operation of 
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a non-trivial georeferencing function to convert back and forth between coordinate 
systems. 

Most FEMA maps do not contain any latitude/longitude information that would be 
useful in georeferencing their scanned raster image. Some recently FEMA-issued panels 
5 do show latitude/longitude values at certain points thereon, but such panels are in a 
distinct minority when considered in the context of the 100,000+ panels that FEMA 
currently provides. 



SUMMARY OF THE INVENTION 

10 Pursuant to the present invention there is provided a system and method for 

georeferencing a digital raster map by associating specific points on the raster map with 
corresponding points on a georeferenced vector map or another georeferenced raster 
map. According to the preferred embodiment, a raster map and a corresponding vector 
map are simultaneously displayed to a user. The user then locates a common 

15 geographic point or feature on each map, and marks each of them as a unique point- 
pair. When the user has marked at least two point-pairs, the system then computes a 
georeferencing function, based on the pixel-coordinates of the points marked on the 
raster map and the corresponding geographic coordinates of the points on the vector 
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map. Thereafter the geographic coordinates of any point on the map may be 
computed. As more point-pairs are identified, the georeferencing function is modified 
for increased accuracy. 

A digital raster map is a computerized map image that resembles a conventional 
5 paper map in that it presents an image of a mapped area, but is not associated with 
any additional underlying data associated with the features of the map. A digital raster 
map is typically created by scanning a conventional paper map, and is a grid-based map 
composed of pixels (or points) of color, shades of grey or black and white. Each pixel 
in the grid has can be referenced by X,Y pixel coordinates, and has only a value, which 

10 indicates the color of that pixel. The pixels of a raster map are not, without more, 
georeferenced, i.e., the pixels of a raster map cannot be identified by latitude/longitude 
or other earth reference. Raster images are commonly referred to as "bit mapped." 

A vector map comprises stored data defining nodes, lines and polygons, rather 
than pixels, to describe an image. Unlike a raster map, in which the map image is 

15 simply stored as a grid of pixels, when a vector map is displayed, it is "drawn" from a 
set of underlying data. A vector map is created through the placements of nodes on a 
plane and connecting those nodes with lines to form closed polygons. Vector lines can 
be attributed with tables of data such as elevations, values, names or other information 
relative to the line. Vector data can be displayed in three dimensions if the lines are 

12 

108344-00010 




attributed with Z values, modified or changed relative to user need, or layered to allow 
for turning off and on the viewing of different information. 

Because of their feature attribution properties, vector maps are particularly 
useful for displaying geographic data. Vector maps are used to display boundaries or 
5 lines that denote the position and extent of features, such as county boundaries or lines 
denoting stream and river systems.^It is.also very easy to view or manipulate the data 
underlying a vector map, for example to view or change the elevation of a feature. 

Also because vector maps are commonly drawn from a geographic data set 
describing the area shown, they are typically georeferenced. Georeferencing is the 
10 process of relating source coordinates, such as X,Y Cartesian coordinates, to geographic 
coordinates, typically in standard latitude/longitude. An image or a vector file is 
georeferenced to be used within a mapping/geographic environment. In a vector map, 
the data from which the map is drawn typically includes a geographic coordinate set. 
Modern GIS systems normally make use of digital vector map information. 
15 However, a vast legacy of paper-based map information exists. It is very expensive and 
time consuming to convert all of the information on these paper maps over to a digital 
vector format. In many cases the scope and expense of such conversions render them 
completely impractical. However, even when a complete conversion to digital vector- 
based format is not possible, it is still possible to obtain some of the benefits of 
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computerized map systems in acordance with the present invention pursuant to which 
paper maps are converted to digital raster maps (by scanning them) and the raster 
image is thereafter georeferenced. After georeferencing, there is, according to the 
present invention, a definite mathematical relationship between the pixel coordinates in 
the raster map, and the geographic coordinates of the feature represented by that 
pixel. The present invention contemplates a system and method for georeferencing a 
raster map by associating points thereon with corresponding points on a previously 
georeferenced vector map. The Invention is described as applicable to georeferencing 
FEMA flood maps, but may be used to georeference any non-georeferenced paper or 
other map. 

First, using known scanning hardware and software, FEMA or other paper map 
images are scanned to create digital raster images thereof. The data representing the 
X,Y coordinates of the pixels comprising the raster images are stored in any convenient 
storage facility, such as computer memory or a CDRom. The stored raster images may 
be referred to herein as "flood map images." 

Second, polygons that approximate the neatline borders of each panel map and 
polygons that approximate the inset boundaries of the flood map images within each 
panel map neatline are "captured." Both types of polygons (and other desired 
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polygons, e.g., those representing notable features on the paper map) are captured, by 
activating a device such as a mouse.. 

The capture of any of the foregoing polygons shown on the raster image, Figure 
2, of a scanned paper map is achieved in the following manner. The user identifies the 
5 vertices (points) of the polygon by "clicking" a mouse icon, in order, on a suitable 
- sequence of points on a displayed flood map image or raster, image. The sequence of., ... 
points may be "smooth" or may be angular, but, in any event, represents an 
approximation of the actual polygon. In the case of the neatlines, because such are 
rectangles, the polygon representation thereof and the actual neatline will be similar, if 

10 not the same. Because inset boundaries are often complexly curvilinear, a captured 
polygon representing an inset is likely to be less "smooth" than the actual polygon. Of 
course, smoothness may be increased by increasing the number of points or vertices 
captures in drawing the polygon. The order in which the user clicks on the vertices 
determines the order of traversal of the vertices of the polygon. As the user identifies a 

15 first and then subsequent vertices, a line segment is drawn from the previous vertex to 
the next identified vertex, until the polygon is completed upon subsequently clicking on 
the first vertex. The completed polygon is preferably displayed as overlying the 
displayed flood map image and may be edited, e.g., by adding additional vertices 
thereto to more accurately mimic the shape of the extant polygon. 
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Elementary tests for the correctness of the polygons thus defined and captured 
are preferably performed. These tests include at least a test to see that no edges of a 
polygon cross. Additional tests may include a test to see that no polygon edge is 
shorter than some edge length parameter, and a test to ensure that the edges defined 
5 do form a closed polygon. These tests are applied on an ongoing basis, as the user 
determines additional vertices of the .polygon. Errors may be.Jiighlighted in different 
colors or brought to the attention of the user through pop-up panels on the computer 
system. The user is allowed to edit the polygons to correct errors, or to improve their 
quality. Each polygon definition is stored, and once the map to which it relates is 
10 georeferenced, the polygons —that is, their vertices and lines— will also be 
georeferenced, 

Next, the flood map (digital raster) images are georeferenced. Initially, there is 
obtained or generated a digital vector map covering the same general geographic area 
covered by the flood map images and capable of displaying and depicting the same 
15 geographic features, such as streets, rivers, railroads, etc. The digital vector map may 
be a commercially obtainable map covering a limited area or an entire geographic 
entity; vector maps covering the entire United States are available. 

Then, the vector map image and the flood map image are simultaneously 
displayed. Such may result in the maps being displayed side-by-side, over-and-under, 
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or in superimposition. Using location information from the flood map image (e.g. 
county name, state, community name, large metropolitan area, etc.), the user 
manipulates the vector map (by zooming, panning, scrolling etc.) until it displays 
roughly the same region displayed in the flood map image. 
5 The user then identifies a first point-pair, that is, a pair of points --one point on 

the vector map, the. other point on a. selected, (if .more than one) inset of the flood 
image map— that represent the same geographic location. Such point-pairs may be 
identified by "clicking" a mouse icon on an identifiable point on the flood map image — 
say, a street intersection— and then "clicking" on the same point or street intersection 

10 shown on the vector map. Accurate placement of the two points can be assured by 
zooming in on the images. Additional point-pairs (one point from the vector map, and 
one from the flood map image) are then identified. 

As the second and subsequent point-pairs are each identified, an algorithm 
calculates and provides a georeferencing function correlating the flood map image and 

15 the raster image. Specifically, upon identification of two point-pairs, a georeferencing 
function is generated and becomes available for the flood map image. This 
georeferencing function is refined as additional point-pairs are identified. As soon as a 
georeferencing function is available it is used to "synchronize" the vector map with the 
digital flood map image. The synchronized maps display the same geographic region, 
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and if one image is scrolled, panned, or zoomed by the user, the same manipulation 
occurs automatically in the other map image. The afore-noted "synchronization" greatly 
facilitates the placement of additional point-pairs, which improves the quality of the 
current georeferencing function, and is more specifically described in commonly 
5 assigned US patent application Serial Number 537,161. 

If multiple inset polygons have been defined on a single flood map image, and if 
those insets are riot shown in their natural correct geographic relationship to one 
another, then it is necessary to calculate a different georeferencing function for each of 
the inset polygons. This situation occurs when FEMA shows a geographic region on one 

10 map panel, which does not actually belong in the same geographic area as the rest of 
the map panel contents. FEMA sometimes does this in order to eliminate the necessity 
of printing another map with only a tiny region shown on it. In contemplation of this 
event, this step of the present invention is more generally stated as follows: When 
georeferencing an inset, the user must select only point-pairs contained within that 

15 inset; a georeferencing function is calculated for each inset, and must be stored so that 
this correspondence is clear. 

Using the foregoing approach, all FEMA flood map insets are georeferenced and 
stored. The georeferencing of the FEMA maps also georeferences the previously 
captured polygons representing the neatlines of the panel maps and the insets 
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therewithin. Accordingly, for every inset on a flood map, there exists a function which 
allows conversion of its X,Y pixel coordinates into lat/lon geographic coordinates. 
Among other things, such georeferencing implies that the boundaries of an inset itself - 
-that is the stored polygon mimicking the inset boundaries™ can be expressed in terms 
5 of lat/lon coordinates. More generally, any identifiable points or regions of a flood map 
image, including the borders of flood zone regions therein, can be expressed in lat/lon 
coordinates. 

If an index map is initially georeferenced, as set forth above> then a preliminary 
georeferencing of a specific regular panel, the neatline of which is shown on the index 

10 map, can be obtained by creating point-pairs out of the corresponding four corners of 
the neatline boundary which is, of course, shown on both the index map and the 
regular panel map. In other words, it is not really required for the georeferencing 
process that one of the maps be a vector map (although this is usually most convenient 
and practical). It is possible for both maps to be scanned raster map images, as long 

15 as the same geographic region is depicted on each map and one map is already 
georeferenced. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The novel features of the invention are set forth in the appended claims. The 
invention itself however, as well as its mode of use, further objects and advantages 
thereof, will best be understood by reference to the following detailed description read 
in conjunction with the accompanying drawings, wherein: 

Figures 1-4 are general representations of four FEMA panel 'maps^thelhsets of 
three of which do not fill the neatlines thereof and the inset of one of which does fill its 
neatlines; 

Figure 5 is a general representation of a FEMA index map made up of a number 
of the panel maps of Figures 1-4; 

Figure 6 is a general representation of a vector map which corresponds to the 
raster map of Figure 2, in accordance with a preferred embodiment of the present 
invention; 

Figure 7 is a flowchart of the process in accordance with the present invention: 

and 

Rgure 8 is a depiction of a computer system for carrying out the method of the 
present invention. 



20 



108344-00010 



DETAILED DESCRIPTION OF THE INVENTION 

The present invention is a system and method for georeferencing digital raster 
maps, Figures 1-5, by associating specific points on a raster map with corresponding 
points on a georeferenced vector map, Figure 6, or another georeferenced raster map. 

Before georeferencing is initiated, the neatlines 32 of each raster map, such as 
Figure 2, are captured as polygons, and the polygons are then stored. Moreover, 
polygons mimicking the border of the inset(s), e.g., 33, 46, 48 in Figures 1, 3 and 4, 
are also captured and stored. The vertices and interconnecting lines of the polygons, 
and preliminary testing methods therefor are discussed above. 

Next, a raster map, such as that shown in Figure 2, and a corresponding vector 
map, Figure 6, are simultaneously displayed to a user. The user then locates a common 
geographic point or feature on each map, Figures 2 and 6. Some corresponding points 
are numbered 200, 202, 204, 206 and 208 on both maps of Rgures 2 and 6. Each set 
of common points 200-200, 202-202, 204-204, 206-206, and 208-208 is "marked" by 
the user as a unique point-pair. Marking may be achieved in any convenient way, 
including "clicking" a mouse on each point of a point-pair so that the x,y coordinates of 
the point on the raster map of Figure 2 and the latitude/longitude coordinates of the 
point on the vector map of Figure 6 are captured and recorded. 

21 

108344-00010 




When the user has marked at least two point-pairs, the system computes a 
georeferencing function, based on the x,y pixel coordinates of the points marked on the 
raster map of Figure 2 and the corresponding geographic coordinates of the points on 
the vector map of Figure 6. Thereafter the geographic coordinates of any point on the 
5 raster map of Figure 2 map may be easily computed. As more point-pairs are 
identified, the georeferencing function is modified for increased accuracy. 

A digital map image is considered georeferenced if a pair of mathematical 
functions, f and g have been defined that can be used to convert back and forth 
between the coordinates of the raster map image of Figure 2 (as defined by the pixels 
10 of the image) and the corresponding longitude and latitude of the location of that point 
on the vector map of Figure 6. That is, georeferencing exists if the following are true: 

1. If (x, y) represents a location on the digital raster 
map image (Figure 2), then f (x, y) =(Lon, Lat) represents 
the longitude and latitude of the corresponding physical 

15 location as depicted on the vector map (Figure 6). 

2. If (Lon, Lat) represents a physical location that lies 
within the region covered by the vector map of Figure 6, 
then g (Lon, Lat) = (x, y) represents the point on the digital 
raster map image of Figure 2 that corresponds to that 

20 longitude and latitude. 
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Here, x and y represent the natural internal coordinate system of the raster map 
image. Typically, as described above, the pixels of a digital raster map, as in Figure 2, 
constitute a natural internal coordinate system. A vector map, as in Figure 6, is based 
on a longitude and latitude coordinate system as its underlying coordinate system. As 
such, a vector map may be considered to be trivially georeferenced. Therefore, the 
functions f and g above .are non-trivial georeferencing functions required to convert 
back and forth between the coordinate systems of the maps of Figure 2 arid Figure 6. 

Referring to Figure 7, first, the raster map, Figure 2, to be georeferenced is 
displayed on a computer monitor at the same time as a previously georeferenced map, 
Figure 6. (step 400). Typically, the map of Figure 6 is a vector map, but could also be 
another raster map for which georeferencing functions are already known and can be 
accessed. The maps, Figures 2 and 6, are preferably displayed side by side, but may be 
displayed above and below, in superposition or in any other simultaneous or pseudo- 
simultaneous (e.g., back-and-forth toggling) display mode. 

The maps of Figures 2 and 6 are preferebly initially independently manipulate 
for rotation, position, pan, zoom and scale. The user is free to manipulate the display of 
both maps until the displayed region of the map in Figure 2 is generally shown in the 
display of the map of Figure 6. 
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Next, pairs of points, or point-pairs, 200-200, 202-202, etc, representing identical 
geographic locations on the maps of Figures 2 and 6 are identified, and the user uses a 
keyboard, mouse, or other pointing or capturing device to mark these point-pairs (step 
405). For each point-pair, 200-200, etc., one point is marked on Figure 2 and the 
5 other point is marked at the corresponding location on Figure 6. The point marked on 
the raster map, Figure 2, has specific x,y coordinates. The same point marked on the - 
vector map, Figure 6, has specific latitude,longitude coordinates. Examples of common 
georeferencing point-pairs include street intersections, 200-200, 202^202, 204-204 and 
206-206, sites 208-208 where a road crosses a river, mountain peaks, specific building 
10 locations, the intersection of two rivers, and other readily identifiable map features. 

If any points on the map of Figure 2 are identified on the map by longitude and 
latitude --in the case of FEMA flood maps, this is a rare occurrence-- a georeferencing 
point-pair can be obtained by marking that point on Figure 2 and marking on Figure 6 
at the location having the identical latitude and longitude. Alternatively, once such 
15 latitude,longitude-identified points on Figure 2 are marked, each will be associated with 
both lat,lon and x,y coordinates, and each can thereafter be used to determine the 
georeferencing functions without requiring a corresponding mark on Figure 6 

After a point-pair is marked on Figures 2 and 6, the x,y coordinate on Figure 2 is 
assigned the corresponding latitude,longitude of its matching point on Figure 6 (step 
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412). When two georeferencing point-pairs are marked, the system computes an 
initial georeferencing function for Figure 2 based on a linear transformation that allows 
an arbitrary rotation and predefined scaling between the x,y domain and the 
latitude,longitude domain (step 415). 
5 To improve accuracy and allow increased generality, more georeferencing point- 

pairs will frequently be desired. To facilitate the process of gathering Jthese additional 
georeferencing point-pairs, the maps of Figures 2 and 6 are synchronized. Accordingly, 
as either map, Figure 2 or Figure 6, is zoomed, panned, scrolled, or otherwise caused to 
display a different region, the other map, Figure 6 or 2, automatically does the same. 
10 This greatly increases the ease with which georeferencing point-pairs can be identified 
and marked on the maps, because, when the user of the system locates an identifiable 
feature on one map, the other map will automatically be showing substantially the same 
location. 

After synchronization is achieved, when the user marks a georeferencing point 
15 on one of the maps, the-system automatically marks the corresponding point on the 
other map in accordance with the last-calculated georeferencing function (step 420). 
The user can then make corrections to the location of a point marked by the system in 
this way (step 425). The user may enter as many point-pairs as are desired; each 
additional point-pair increases the accuracy of the georeferencing process. 
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When three or more georeferencing point-pairs are determined, the system 
computes a georeferencing function for Figure 2 based on a general linear 
transformation between the x,y domain and the latitude,longitude domain. When four 
or more georeferencing point-pairs are determined, the general linear georeferencing 

5 function is over-determined. This means that more than the required amount of 
information to compute the general linear georeferencing functions is available, but that 
this information may not be completely consistent. The system uses the extra 
information contained in the additional georeferencing points to provide validation 
checks to protect against the possibility that some of the data points may be inaccurate 

10 (step 430). Point-pairs that deviate excessively with respect to a calculated standard 
error are presumed to be inaccurate and are omitted from the extant calculation of the 
georeferencing function. As new point-pairs are added, the system rechecks previously 
omitted points to determine if those previously omitted points should now be 
considered when recomputing the georeferencing functions. 

15 The additional point-pairs allow the system to compute the general linear 

georeferencing functions which best fit the combination of all the available 
georeferencinOg point-pairs, so that the active georeferencing functions are revised 
with each new point-pair (step 435). According to preferred embodiments, this is done 
by using a "least square" parameter fitting operation. 
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The user may then proceed to enter the next point-pair (step 440) . When the 
user is finished, the system stores the active georeferencing functions with the raster 
map of Figure 2 (step 445). At this time, the raster map is considered fully 
georeferenced and the related georeferencing functions are stored as is their 
5 association with the Figure 2. When the raster map is later accessed, the system 
retrieves the relevant georeferencing functions, and applies them to determine the 
latitude and loOngitude of any point on the raster map or to determine the x,y 
cordinates of any point on the vector map of Figure 6. Thus, at this time, the neatlines 
of each panel map, the boundary of each panel map inset and the "grid" of neatlines on 
10 each index map are, as polygons, fully georeferenced. 

While various techniques exist for determining a georeferencing function from 
sets of point-pairs, the specific approach used by the system and method of the present 
invention is now discussed. 

To carry out the georeferencing process of the present invention, sets of points 
15 P, are specified. Let / e P denote one of the points i of P . Associated value x if y t and 

Lat,,Lon t refer, "respectively, to the* and y coordinates of the point i on the raster 

map of Figure 2, and the latitude and longitude of the same point i on the earth's 
surface, i.e., on the vector map of Figure 6. There is also defined a set of points Ac P, 
where A is the subset of "active" points, i.e. points which are in extant use in the 
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calculation of the georeferendng functions. A parametric family of functions □ is 
selected. It is stipulated that functions in this family □ map from (x,y) to(Lon,Lat). 

Any function / e □ has components f l (x,y)=Lon, and f 2 (x,y) = Lat. From within this 
family there is sought the specific function /which (in a sense to be defined below) 
5 introduces the smOallest error in attempting to render the following system of equations 
true: 

f(x i ,y) = (Lon i ,Lat) for is A (1) 
Once determined, /will function as the georeferendng function which is used to 
compute corresponding latitude and longitude values (Lon,Lat) on the vector map of 

10 Figure 6 for any point (x,y) on the raster map of Figure 2. There are a variety of ways 
to define the function that minimizes error in rendering Equation (1) true. The method 
of the present invention utilizes a least squares approach, known in mathematics as an 
L 2 approach. This approach seeks to ascertain the function / which minimizes the sum 
of the squared differences between the actual and the predicted values of latitude and 

15 longitude. In other words, from among all the functions/ e □, /is the one which 
minimizes: 

SSE = jfoy-W +{f 2 (x i ,y l )-Lat t y\ (2) 
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Alternatively, the function /may be chosen so that the sum of absolute errors, 
rather than squared errors, is minimized or so that the largest error is minimized. Other 
criteria are also possible. 

In what follows there shall be described methods for computing georeferencing 
functions based on both a general linear parametric family of functions and a rotational 
linear family of functions. Those skilled in the art will readily perceive that this approach 
may be easily generalized to other parametric families of functions. One such other 
parametric family is a nonlinear parametric family of functions arising from the 
projection of the earth's surface onto a flat map. 

The General Linear Case 

In the general linear case, F is the set of all possible linear transformations 
which map from (x,y) to(Lon,Lat). Thus, 

*1 (3) 

Pi. 

for some choice of parameters 

/ ^12/ ^21/ ^22' ^i/ b 2 ' If the region 

covered by the map to be georeferenced is not too large, then this family of functions 

will contain a suitable function / whose total error, as determined by the least squares 
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methodology, is quite small. In the case where the map to be georeferenced covers an 
area that is too large, the curvature of the earth must be taken into account and F is 
not a suitable family of functions. In the latter event, nonlinear functions must be used 
as mentioned above. 

To find / we seek the parameters which minimize 



SSE = £ (a n x, + a u y, +b t - Lon] f + (d 2l x, + a^y, + b 2 - Laf,) 2 



(4) 



ieA. 



The parameter values which minimize this expression are found by solving the 
following two independent systems of linear equations: 
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These systems can be solved by well known methods, such as Gaussian 
Elimination or LU factorization. The solutions yield the desired values ofa„, d l2f <5 21 , 

a al b\, and b 2 . It should be noted that equations (5a) and (5b) do not have a unique 
solution unless three or more non-colinear points on the raster map, Figure 2, are 
contained in A . Generally speaking, then, it requires the selection of three point-pairs 
to initiate the calculation of a georeferencing function from the family of general linear 
transformations. When there are four or more point-pairs, it is possible to compute a 
standard deviation of errors using the formula: 



where s is an estimator for the amount of error to be expected between actual 
and predicted latitude and longitude values. 

The inverse georeferencing function, /"' that maps from (Lon, Lat) of one point 
of a point-pair on the vector map of Figure 6 to (x,y) of the other point of the point-pair 




(6) 
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on the raster map, Figure 2, may now be readily obtained by inverting the function / . 
Having done this, it is possible to compute a similar standard deviation of error for 
which is an estimate of error as measured in the x,y coordinates of the raster map of 
Figure 2. 

The Rotational Linear Case: 



10 



In the rotational linear case, F is the set of all possible linear transformations 
which map from (x,y) to (Lon,Lat) and which also allow for a known scale change 
caused by the difference in distance per degree of latitude and per degree of longitude 
at any point not on the equator. Thus 
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(7) 



for some choice of parameters p x , p 2l p if and p AI where y =cos(Lat) is 
computed from the estimated latitude of the map to be georeferenced. Clearly, this is a 
special case of the more general linear transform described above, where 
15 a„ = P s / r,a 12 = -p 4 ,a n = -p 4 , 
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and 

In theory, if the region covered by the raster map of Figure 2 to be 
georeferenced is not too large, then this family of functions will contain a suitable 
5 function, / whose total error is quite small. In the case where the map to be 
"georeferenced covers a larger area than this, the curvature of the earth must be taken 
into account and F \s not a suitable family of functions. In such a case, nonlinear 
functions must be used, as noted earlier. In practice, the general linear transform will 
yield results which are somewhat superior to the rotational linear transform due to the 
10 former's innate ability to use arbitrary scaling factors. These arbitrary scaling factors 
might compensate, for example, for scaling flaws in the digital scanner used to scan in 
the non-georeferenced paper maps to produce maps like that of Figure 2. The primary 
advantage of the rotational linear transform is that it can be computed with only two 
reference point-pairs, while the general linear transform requires three reference point- 
15 pairs. 

To find /there is sought, as before, the parameters which minimize .SSEas 
defined in equation (4). 
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The parameter values which minimize this expression are found by solving the 
following system of linear equations: 

« ° z*/ -z^.- ] r fz, L ° n i 

leA /eA leA „ n /? ,€ ^U 

. /eA /eA /eA J L <eA /eA ... 



These systems can easily solved by well known methods, such as Gaussian 
Elimination, or LU factorization. The solutions yield the desire values of y9, , fi u p x , 
and fi x which in turn yield the desired values for <5„, d u , d 21l d a , b t , and b 2 . 

It should be noted that equation (8) does not have a unique solution unless two 
or more points are contained in A. Generally speaking, it requires two points to 
determine a georeferencing function from the family of rotational linear 
transformations. When there are three points or more, it is possible to compute a 
standard deviation of error, s using the formula: 
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s = 




(9) 



where s is an estimator for the amount of error to be expected between actual 



5 



and predicted latitude and longitude values. 



The inverse georeferencing function / _1 that maps from (Lon, Lat) to (x,y) is 



a similar standard deviation of error for which is an estimate of error as measured 
in the bitmap coordinates. 



When point-pairs are selected, the apparent values of x,y and Lon,Lat which are 
selected by mouse "clicking" may include error. To reduce the risk of incorrect 
georeferencing resulting from such errors, certain error handling procedures are built 
into the georeferencing process of the present invention. One technique utilized by the 
15 present invention is to detect a "bad" point-pair and then remove it from the set A of 
"active" point-pairs, that is, those point-pairs presently involved in the calculation of a 
georeferencing function. Removing a bad point-pair from A does not delete the 
information associated with that point-pair. Rather, removing the bad point-pair causes 



readily obtained by inverting the function / . Having done this, it is possible to compute 



10 



Automatic Error Detection and Handling 
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the extant calculation of the georeferencing parameters to be uninfluenced by that 
point-pair. TTie point-pair is not removed entirely, because it may be determined at a 
later stage of the georeferencing, that the point-pair was not really bad, and should be- 
used in the georeferencing calculation. This will be clarified shortly. 

The following steps outline the bad point detection process using the general 
linear transform approach to georeferencing. 

(1) Place all point-pairs so far selected into the active set A. An "active" 
point-pair comprises a pair of points which are presently being used in the extant 
georeferencing calculation. 

(2) If there are four or less active point-pairs, error correction cannot be 
performed and the georeferencing function cannot be refined. 

(3) If there are five or more active point-pairs, temporarily exclude each of 
them in turn from the active set A, designating a temporarily excluded point-pair as k. 
With each point k excluded in turn, the calculate the inverse georeferencing function 
£ (t) using the non-k point pairs. Calculate SSE k for g 



(4) Calculate the difference between the predicted value and the actual value 
6 k -]p*\Lon t ,Lat k )-(x kt y k ]. 



(5) Save the values, S k and 
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(6) Return point k to the active set (and move on to the next value of k . 

(7) From among the results found in steps (3)-(6) / ascertain if there exists a 
point k , the absence of which from the calculation of the inverse function, leads to the 



largest values of y SSE >c, and S k , where the value of 5 k also satisfies 
5k /Co F )c t andS k )c 2 , where c,, andc 2 are constants which are set according to the 



general level of accuracy to be expected on the particular type of map which is being 
georeferenced, the current number of active point-pairs, and the dots or pixels per inch 
of the raster image. If there is such a point, it is designated as erroneous and is 
removed from the active set A, and return to steps ( 2) and (3) above. If there is no 
such point, error detection is completed. 

The foregoing error detection may be, of course, carried out on a general 
purpose computer, preferably the same computer included in the system for 
manipulating the various maps and for calculating the georeferencing functions. 

There are several things to note about the foregoing error detection procedure. 
One is that allowing the value of c, and c 2 to change with the number of active point- 
pairs, makes it possible for the georeferencing system and method to utilize point-pairs 
which might originally have been determined as erroneous or inconsistent after a large 
enough sample of points has been gathered to make it clear that a lesser level of 
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accuracy is all that can be achieved on this map. Use of the error detection procedure 
of the present invention cannot reduce the number of active point-pairs to less than 
four. This scheme means that as each new point is added, all points determined so far 
are considered, even those which had previously been marked bad. Thus early 
5 "misjudgements" on the part of the system can be corrected later, in light of new point 
information. 

The same bad point detection process, can also be implemented when the 
rotational linear transform is used. In this case the method is capable of reducing the 
number of active point-pairs to a minimum of three, rather than four for the general 
. 10 linear transform approach outlined above. This can be useful when dealing with small 
sets of active points. 

When the system is georeferencing raster map images that cover a large enough 
area so that the (nonlinear) curvature of the earth is a source of significant error, the 
system can calculate nonlinear georeferencing functions whose form corresponds to the 
15 map projection that was used to create the map of Figure 2. 

A specific example of the operation and application of the preferred 
georeferencing method may be shown with reference to the "Flood Zone 
Determination" business. The Federal Emergency Management Agency (FEMA) 
publishes a library of tens of thousands of paper maps showing various types of flood 
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zones and their locations in the United States. A flood zone determination on a property 
is frequently done in the following way: 

1. The address of the property is examined, and the 
location of the property is determined (perhaps through the use of 
a geocoding system, or by examining an available street map). 

2. A map analyst attempts to determine which of the 
many thousands of FEMA flood maps will contain this property. 

3. The map analyst goes to a map storage area and 
retrieves the desired map, often examining several maps before 
making a final selection. 

4. Having retrieved the paper map, the map analyst 
next determines where, precisely, the property is located on the 
map. 

5. Finally, the map analyst examines flood zone 
notations on the map at the property's location in order to 
determine its flood-zone status. 
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When performed using paper maps, the above process is difficult and quite time 
consuming . A database of scanned raster map images, alone, can be used to reduce 
the time and effort expended on step 3 above. However, georeferenced raster map 
images can also be used to eliminate virtually all work from steps 2 and 4 above. Since 
5 those are generally the most time consuming steps under the current process, the value 
of georeferenced flood maps is considerable. 

Using georeferenced flood map raster images, steps 2 and 4 above, are replaced 

by: 

2'. A computer system combines the pre-designated outlines of the raster map 
10 and the georeferencing information to obtain a polygon expressed in terms of latitude 
and longitude that outlines the region included in each flood map. Then the system 
determines which of the polygons contain the address in question, which is done 
using a w point-in-polygon" algorithm. At the conclusion of this process, the computer 
system has identified a map panel (or perhaps a small number of map panels) that 
15 contains the address. 

4'. Since the latitude and longitude of the property are known (by virtue of a 
geocoding phase), the computer system can use the georeferencing of the map 
panels to locate the property on each of the panels found above, thus largely 
eliminating any need for the map analyst to scan the flood map for the address location 
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Therefore, it is easily seen that by using georeferenced raster maps instead of 
paper maps, the process of determining flood zones is substantially automated and 
much more efficient. 

It is important to note that while the present invention has been described in the 
5 context of a fully functional data processing system and/or network, those skilled in the 
art will appreciate that the mechanism of the present invention is capable of being 
distributed in the form of a computer usable medium of instructions in a variety of 
forms, and that the present invention applies equally regardless of the particular type of 
signal bearing medium used to actually carry out the distribution. Examples of 
10 computer usable mediums include: nonvolatile, hard-coded type mediums such as read 
only memories (ROMs) or erasable, electrically programmable read only memories 
(EEPROMs), recordable type mediums such as floppy disks, hard disk drives and CD- 
ROMs, an transmission type mediums such as digital and analog communication links. 
While the invention has been particularly shown and described with reference to 
15 a preferred embodiment, it will be understood by those skilled in the art that various 
changes in form and detail may be made therein without departing from the spirit and 
scope of the invention. 
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WHAT IS CLAIMED IS: 

1. A computer-implemented method of georeferencing an initially non- 
georeferenced raster map image, comprising the steps of: 

(a) displaying a raster map and a georeferenced map that both depict the 
5 same geographic area; 

(b) identifying at least two point pairs, each point pair comprising one 
point on the raster map and one point on the georeferenced map, each point of a point 
pair lying on the same geographic feature on both maps; 

(c) with the use of the raster image coordinates of each point on the 
10 raster map and the geographic coordinates of the corresponding point on the 

goereferenced map, determining the functional relationships between the image 
coordinates and the geographic coordinates; 

(d) thereafter, for each additional corresponding point-pair that is 

identified 

15 (i) revising the functional relationships between the image 

coordinates and the geographic coordinates according to the additional point pairs, and 

(ii) disregarding any points which are substantially inconsistent 
with the functional relationships. 

2. A method as in Claim 1, wherein: 

20 the image coordinates are Cartesian (x,y) coordinates, and the 

georeferenced map is a vector map having (latitude,longitude) coordinates. 

3. A method as in Claim 1, wherein: 

the maps are displayed substantially simultaneously in step (a). 
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4. A method as in Claim 1, wherein: 
although the maps diversely depict the geographic area, certain 
geographic features and their intersections are displayed in common thereon. 



5 5. A method as in Claim 4, wherein: 

the identical geographic features include roads, streets, stream beds, 
railroad tracks, mountain peaks and buildings. 

6. A method as in Claim 1, wherein: 

in step (c), the functional relationships between the image coordinates 
10 and the geographic coordinates of the point pairs are calculated by selecting from a 
family of parametric equations the one equation that introduces the smallest error 
between the actual and predicted values of longitude and latitude assigned to (x,y) in 
the image. 

7. A method as in Claim 6, wherein: 

15 the one equation that introduces the smallest error is determined 

according to a least squares technique. 

8. A method as in Claim 6, wherein: 

the selection of an equation from a family of parametric equations is 
effected according to a general linear transform. 
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9. A method as in Claim 8, wherein: 

at least three non-colinear point-pairs are Identified. 

10. A method as in Claim 6, wherein: 

the selection of an equation from a family of parametric equations is 
effected according to a rotational linear transform. 

11. A method as in Claim 10, wherein: 

at least two non-collinear point-pairs are identified. 

12 A method as in Claim 1, further comprising: 
storing the functional relationship with the raster map. 

13. A method as in Claim 1, wherein: 

the point pairs are identified in step (b) by a user marking them following 
a visual determination that the geographic features on which they lie are the same 

14. An item of commerce, comprising: 

a computer-usable, portable memory unit containing the data of Claim 12. 

15. An item of commerce, comprising: 

a data transmission medium carrying the data of Claim 12. 
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16. An item of commerce, comprising: 

a memory unit containing computer instructions for effecting the method 

of Claim 1. 

17. An item of commerce, commerce: 

5 a data transmission medium connected to a general purpose computer for 

instructing the computer to effect the method of Claim 1. 



18. A computer system for georeferencing an initially non-georeferenced 
raster map image, comprising: 

10 (a) a monitor for displaying a raster map and a georeferenced map that 

both depict the same geographic area; 

(b) first facilities for identifying at least two point pairs, each point pair 
comprising one point on the raster map and one point on the georeferenced map, each 
point of a point pair lying on the same geographic feature on both maps; 

15 (c) second facilities that use of the raster image coordinates of each point 

on the raster map and the geographic coordinates of the corresponding point on the 
goereferenced map, for determining the functional relationships between the image 
coordinates and the geographic coordinates; 

(d) third facilities for, after each additional corresponding point-pair is 

20 identified, revising the functional relationships between the image coordinates and the 
geographic coordinates according to the additional point pairs, while disregarding any 
points which are substantially inconsistent with the functional relationships. 
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19. A system as in Claim 18, wherein: 

the image coordinates are Cartesian (x,y) coordinates, and the 
georeferenced map is a vector map having (latitude,longitude) coordinates. 

20. A system as in Claim 18, wherein: 

5 the maps are displayed substantially simultaneously by the monitor. 

21 A system as in Claim 18, wherein: 

although the maps diversely depict the geographic area, certain " 
geographic features and their intersections are displayed in common thereon. 

22 A system as in Claim 21, wherein: 

10 the identical geographic features include roads, streets, stream beds, 

railroad tracks, mountain peaks and buildings. 

23. A system as in Claim 18, wherein: 

the second facilities calculate the functional relationships between the 
image coordinates and the geographic coordinates of the point pairs by selecting from a 
15 family of parametric equations the one equation that introduces the smallest error 
between the actual and predicted values of longitude and latitude assigned to (x,y) in 
the image. 

24. A system as in Claim 23, wherein: 

the one equation that introduces the smallest error is determined 
20 according to a least squares technique. 
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25. A system as in Claim 23, wherein: 

the selection of an equation from a family of parametric equations is 
effected according to a general linear transform. 

26. A system as in Claim 25, wherein: 

5 at least three non-colinear point-pairs are identified. 

27. A system as in Claim 23 wherein: 

the selection of an equation from a family of parametric equations is 
effected according to a rotational linear transform. 

28. A system as in Claim 27, wherein: 

10 at least two non-colinear point-pairs are identified. 

29 A system as in Claim 18, further comprising: 

fourth facilities for storing the functional relationships together with the 

raster map. 

15 30. A system as in Claim 18, wherein: 

the first facilities include a device by which a user can mark the points of a 
point pair following a visually determining from the monitor that the geographic 
features on which the points lie are the same. 



47 



108344-00010 



# 



ABSTRACT OF THE DISCLOSURE 



After FEMA paper flood maps are scanned into a computer system, they are each 
displayed side-by-side with the relevant portion of a stored vector map covering the 
geographic area covered by all of the FEMA maps. During display, corresponding 
5 points, or "point-pairs," on each map are marked by a keyboard or mouse device. The 
corresponding points may be such artifacts as roads, stream beds, railroad tracks, 
intersections of the foregoing, buildings, mountain peaks, etc., which are identifiable on 
each displayed map. When two point-pairs are marked, a stored algorithm calculates 
georeferencing functions which permit translation of the x,y coordinates of the raster 
10 map to the latitude,longitude coordinates of the vector map and vice versa. As more 
point-pairs are marked, the georeferencing functions are refined. Another algorithm 
removes apparent "bad" points from the georeferencing calculations, which points may 
re later used if it appears that they are not bad. 
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